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ABSTRACT 

We report on the VLBI follow-up observations using the Japanese VLBI Network (JVN) array at 
22 GHz for the largest X-ray flare of TeV blazar Mrk 421 that occurred in mid-February, 2010. The 
total of five epochs of observations were performed at intervals of about 20 days between March 7 
and May 31, 2010. No new-born component associated with the flare was seen directly in the total 
intensity images obtained by our multi-epoch VLBI observations. However, one jet component located 
at ~ 1 mas north-west from the core was able to be identified, and its proper motion can be measured 
as —1.66 ± 0.46 mas yr _1 , which corresponds to an apparent velocity of —3.48 ± 0.97c. Here, this 
negative velocity indicates that the jet component was apparently moving toward the core. As the 
most plausible explanation, we discuss that the apparent negative velocity was possibly caused by the 
ejection of a new component, which could not be resolved with our observations. In this case, the 
obtained Dopplcr factor of the new component is around 10 to 20, which is consistent with the ones 
typically estimated by model fittings of spectral energy distribution for this source. 

Subject headings: BL Lacertae objects: individual (Mrk 421) — techniques: interferometric 



1. INTRODUCTION 

Blazars are a sub-class of radio-loud active galac- 
tic nuclei (AGNs) with relativistic jets pointing almost 
along the line of sight. The spectral energy distribu- 
tions (SEDs) of blazars are dominated by non-thermal, 
strongly Doppler boosted, and variable radiation pro- 
duced at the innermost part of the jets. The SED and 
the multi-frequency correlations of blazars have been in- 
tensively studied in the past through dedicated multi- 
frequency campaign observations (eg. Abdo et al. 2011, 
and references therein) 

One of the most important pieces for understanding 
relativistic jet formation is to obtain the direct image of 
a region radiating high energy emission. However, it is 
difficult to obtain these images for the following reasons. 
(1) Many of the radio observations have been done by 
single-dish telescopes (eg. Raiteri et al. 2009; Villata et 
al. 2009) and their spatial resolution is much larger than 
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the relevant size. (2) Recent progress of VLBI obser- 
vations indicate possible correlations between newborn 
VLBI components and high energy emissions by the MO- 
JAVE (Monitoring Of Jets in Active galaxies with VLBI 
Experiment) project (eg. Lister et al. 2009a; Kovalev 
et al. 2009). It is, however, difficult to confirm the 
correlation between a newborn VLBI component and X- 
ray and/or the very high energy (VHE) 7-ray flares when 
several flares occur during a single monitoring span which 
is typically a few months for each source. Therefore, to 
confirm the fact mentioned above, to examine the behav- 
ior at the root of the jet by conducting a dense monitor 
using VLBI soon after (< a few weeks) the giant X-ray 
and/or VHE 7-ray flare, can be a complementary ap- 
proach. 

Mrk 421 is one of the best sources for studying the 
most compact regions in blazars, because of its proximity 
(z = 0.031, 131 Mpc). With respect to the kinematics 
of the parsec-scale jets in Mrk 421, very slow and/or 
almost stationary components, which can be monitored 
by long-term observations with VLBI, have been studied 
by Piner et al. (2010). However, no dense monitoring 
has yet been conducted to study the behavior of near the 
core of Mrk 421 right after the giant flares (see, §2). 

Throughout this paper, we adopt the cosmological pa- 
rameters of Hq = 71 km s~ 4 Mpc -1 , O m = 0.27, and 
fl\ = 0.73 (Komatsu et al. 2009). Using these param- 
eters, 1 mas corresponds to a linear distance of 0.64 pc, 
and a proper motion of 1 mas yr~ 4 corresponds to an 
apparent velocity of 2.1c at the distance of Mrk 421. 

The organization of this paper is as follows: In §2, 
we summarize past giant flares of Mrk 421, and their 
follow-up observations using VLBI. In §3, §4, and §5, we 
describe our observations, data analysis, and the results 
obtained from model-fitting the data. In §6, we compare 
our result with the one derived from the MOJAVE pro- 
gram to identify the position of the jet component, and 



2 



Niinuma ct al. 



in §7, we discuss the astrophysical implications derived 
from our result. 

2. REVIEW OF GIANT VHE 7- RAY AND X-RAY FLARES 

As mentioned in §1, systematic VLBI monitor within 
^100 days for investigating the behavior of near the core 
after the occurrence of a giant flare have not done yet. In 
this section, we briefly review the previous giant flares of 
Mrk 421 at VHE and/or X-ray band, and the follow-up 
observations with VLBI for the flares. 

2.1. Previous flares 

Since the first discovery of VHE 7-ray emission from 
Mrk 421 was reported by Punch et al. (1992), large 
X-ray or VHE 7-ray flares of Mrk 421 have occasion- 
ally been reported. Gaidos et al. (1996) succeeded in 
detecting the VHE 7-ray flares from Mrk 421 within a 
few years after the discovery of VHE emission of Mrk 
421. An extremely high state of Mrk 421 at VHE band, 
which lasted from 2001 January to March was reported 
by Krennrich et al. (2002). Recently, multi-band obser- 
vations have also been performed for such large flares of 
Mrk 421 (BLazejowski et al. 2005; Lichti et al. 2008; 
Donnarumma et al. 2009) 

In Table 1, we listed the large flares of Mrk 421, and 
their VLBI follow-up observations. (Here, we excluded 
the information of geodetic observations). The epoch, 
peak fluxes at X-ray/VHE band, and references of each 
flare are indicated in column 1-4, and the information 
about VLBI observations for these flares is shown in 
column 5-7, if the observations were done within three 
months after the flare occurred. Columns 5 and 6 show 
their first epoch (as the earliest one), and the averaged 
interval of each epoch in the case that multi-epoch obser- 
vations were done, respectively. For example, although 
Piner & Edwards (2005) carried out the multi-epoch 
observations for the exceptionally strong TeV gamma- 
ray flare in 2001 with VLBA, their first epoch was done 
about five months after the flare and the intervals of each 
observation were about three months. 

2.2. The flares detected by MAXI 

The total of four strong flares (> 100 mCrab) of Mrk 
421 were detected with the energy ranging from 2 keV 
to 10 keV between 2009 November 1 and 2010 February 
16 (Isobe et al. 2010) by the Monitor of All-sky X-ray 
Image (MAXI; Matsuoka et al. 2009), on the Exposed 
Facility of the Japanese Experiment Module " Kibo" at- 
tached to the International Space Station (ISS). The sec- 
ond and fourth ones were also clearly detected with the 
energy range of Swift /BAT (15 keV - 50 keV). In the 
fourth one, the flux reached 164 ± 17 mCrab on a 6-hour 
average which is the largest one at X-ray band among 
those ever reported. Especially for the largest X-ray 
flare, the detection of very bright VHE 7-ray emission 
(~ 10 Crab) was also reported (Ong 2010). 

3. OBSERVATIONS AND DATA REDUCTION 

Anticipating the detection of a newly born component 
associated with the largest X-ray flare, we started follow- 
up observations with the Japanese VLBI network (JVN: 
Doi et al. 2006) 11 within a month after 2010 February 

11 http:/ /www. astro. sci.yamaguchi- 

u. ac.jp/~kenta/jvnhp/cng/indcx_c. html 



flare. JVN includes VLBI Exploration of Radio Astrom- 
etry (VERA: Kobayashi ct al. 2003), and is spread 
across the Japanese islands. Our first JVN observation 
of Mrk 421 took place on 2010 March 7. It was followed 
by another 4 epochs distributed by regular intervals of 
~ 20 days. The last observation was conducted on 2010 
May 31. All these observations were performed at 22 
GHz. Unfortunately the data we could analyze are only 
from VERA antennas because the other antennas mainly 
suffered from recording trouble. 

The data from 4 telescopes ("VERA Mizusawa", 
"VERA Iriki", "VERA Ogasawara", and "VERA 
Ishigaki-jima" : Baseline ranging over 1000 — 2300 km) 
that were recorded with a bandwidth of 256 MHz using 
2-bit samples and left-circular polarization, were used 
for data reduction. The on-source time of Mrk 421 in 
each epoch is approximately 5.5-hr, and the range of el- 
evation at which we carried out observations for Mrk 
421 were 35° < El < 88° until the fourth epoch, and 
10° < El < 82° at the fifth epoch. The elevation depen- 
dence of the aperture efficiency of each VERA antenna 
is almost flat at El > 25° 12 . The synthesized beam size 
in our observations is typically a major axis of 1.2 mas 
and a minor axis of 0.8 with a position angle of —50° 
measured from north through east. The correlation was 
carried out on the Mitaka FX correlator (Chikada ct al. 
1991). 

In the last epoch May 31, fringes to the "VERA Mizu- 
sawa" antenna could not be detected due to system trou- 
bles. This means the maximum baseline length was lack- 
ing. Therefore we only had data from three antennas in 
this epoch. Since the synthesized beam size is 1.5 times 
larger than in the other four epochs, we decided to ex- 
clude this observation from further analysis. Also, at 
the second epoch Apr 1, Mrk 421 was not detected in 
the more than half of an observing time at the baseline 
including VERA Iriki station because of bad weather. 

For the correlated data, we performed the standard 
VLBI calibration and fringe-fitting procedure using the 
NRAO Astronomical Image Processing System (AIPS) 
software package. The flux calibration done by the AIPS 
task APCAL typically achieved an accuracy of 10%, for 
VERA at 22 GHz (Nagai et al. 2010). The calibrated 
visibilities were then exported to carry out an imaging 
procedure using the Caltech Difmap package. After the 
time averaging of visibility data in 30-second bins, we 
carried out the iteration of CLEAN and self-calibration 
procedures in Difmap. During this iteration process, 
the solution interval of an amplitude self-calibration was 
shortened from intervals as long as the whole observa- 
tional time down to 15 minutes. 

4. MODELFIT 

The total intensity images of Mrk 421 are shown in 
Figure 1. These images were convolved with a circular 
beam with an axis size (full width at half maximum) 
of 0.8 mas (shown in the bottom left corner of the fig- 
ure), which corresponds to the typical minor axis of a 
synthesized beam in our observations. Our observations 
confirmed a north-west oriented core-jet structure of Mrk 
421. The direction of the outflow is consistent with those 
reported in other studies (Lico et al. 2012; Piner et al. 

12 http: / / veraserver.mtk.nao.ac.jp/ restricted / CFP201 1 / statusl 1 .pdf 
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1999; Pincr & Edwards 2005; Piner et al. 2010). 

Following Pincr et al. (2010), we have carried out the 
elliptical and the circular Gaussian model-fitting to the 
self-calibrated visibility data in the (u — w)-plane using 
the task modelfit in the Difmap at all epochs. In such 
a way as to avoid the unrealistic changes in Gaussian 
sizes between adjacent epochs, we found that the ellip- 
tical Gaussian model for the core and the circular one 
for the jet component named JC1 is an appropriate com- 
bination. Model parameters are given in Table 2. The 
la errors or la upper limits of each parameter were de- 
rived from the model-fit procedure. The parameter we 
focused on is varied around the best-fit value by confirm- 
ing a change in the x 2 value. The position angle of the 
core components is around —30°, which coincides with 
that between the core and JC1 (about —40° except for 
the 2nd epoch). From our model-fit, the model compo- 
nents for the core were resolved toward the jet axis at 
all epochs. Regarding the direction perpendicular to the 
position angle of the core component, one is resolved at 
only the 3rd epoch, and others are unresolved. When 
we fixed the axial ratios of the core components as the 
value obtained at the 3rd epoch, a slight difference be- 
tween the major axis of the core component and the jet 
axis generates a negligibly small fitting error of 0.02 mas 
at maximum. Therefore JC1 is sufficiently distinguished 
from the core at all epochs, because of the size of its 
major axis. 

5. RESULTS 

The position of JC1 relative to the core is plotted as 
a function of time in Figure 2. As the result of the 
least square fit to the displacements of the position of 
JC1 relative to the core during 66 days (from Mar 7 to 
May 12) under the assumption of linear motion of JC1, 
we obtained the proper motion of /ii a = —1.66 ± 0.46 
mas yr -1 which corresponds to the apparent velocity of 
vi a = -3.48±0.97cat the distance of Mrk 421. This neg- 
ative velocity implies apparent inward motion of JC1. To 
measure the apparent velocity of JC1, although we have 
to estimate not only thermal, but systematic errors as 
the uncertainty on the position of JC1 relative to the 
core at each epoch, it is quite difficult to estimate the 
latter quantitatively due to its non-linearity. Therefore, 
following Homan et al. (2001); Pincr et al. (2010), 
we set the positional uncertainty on the data points that 
were uniformly rescaled so that the reduced x 2 1S unity. 
Also, as seen in Figure 2, if we exclude the last epoch 
from the fit, the result obtained would only be half of 
the value mentioned above. Here, we applied the first 
order fit to the motion of JC1. Based on our results, the 
significance of the detected superluminal motion of JC1 
in Mrk 421 is 3.57ct. 

On the other hand, though the 10% flux errors men- 
tioned in §3 are slightly large relative to the core flux 
change, which can be seen in Table 2, it appears that 
the core flux in the last epoch is higher than in the pre- 
vious epochs. Additionally, the other flux data, which 
were compiled by F-GAMMA program 13 (Angclakis et al. 
2010), also shows almost the same trend of increasing 
flux as our result. Although the scale of flux variability 
was not large, there is a possibility that a delayed in- 

13 http:/ /www. mpifr-bonn.mpg.de/div/vlbi/fgamma/rcsults. html 



crease in radio flux was seen. When the flare occurs in 
the blazars, it has often been reported in the radio ob- 
servations at low frequency because of the opacity effect 
(e.g., Orienti et al. 2011). These support the possibility 
of structural changes occurring in the core region, such 
as the emergence of a new component. The details of the 
flux density change at radio band after the X-ray flare of 
2010 February, will be discussed in a forthcoming paper. 

6. IDENTIFICATION OF JC1 

To identify the correspondence between the jet compo- 
nent JC1 and the one detected by other observations, we 
analyzed the data observed with VLB A at 15 GHz (MO- 
JAVE program: Lister et al. 2009c). Wc obtained the 
calibrated visibility data at a total of five epochs (2009 
Dec 17, 2010 Feb 11, Jul 12, Oct 15, and Nov 29) from the 
MOJAVE web site, and carried out the standard analy- 
sis procedure. In Figure 4, we show a CLEANed image 
within 4 mas from the core at the epoch of 2010 Febru- 
ary 12th. Also, model fit parameters are shown in Table 
3. Here, we carried out the model fit in the same way 
as mentioned above. As the result, we identified four jet 
components in this scale by comparing the angular dis- 
tance from the core, the position angle, the size, and the 
flux density of each model throughout five epochs, and 
named them MCI, MC2, MC3, and MC4 from outermost 
to innermost component. In addition, as temporary com- 
ponents, MC3a and MC4a are seen in the third and fifth 
epoch, respectively. 

The angular distance from the core for all model com- 
ponents at each epoch derived from the MOJAVE data, 
and the one detected by our observations (JC1) are 
shown in Figure 5. From the model fit parameters, we 
can identify the JC1 as MC3. MCI, MC2, and MC4 
cannot be detected in the JVN images (e.g., Figure 3). 
The reason MCI and MC2 were not detected is that the 
image sensitivity of our observations was insufficient to 
detect both components which arc relatively faint. And 
MC4 cannot be resolved with the angular resolution of 
our observations either, because of its closeness to the 
core. 

Each solid line in Figure 5 indicates the least square 
fit to the motion of each component using first order fit. 
These fitting were done without JC1 data, and each er- 
ror bar was calculated in the same way as in Figure 2. 
For MC3a and MC4a, each positional error is set the 
same value as MC3 and MC4 respectively, because each 
component was not detected at other epochs. In this fit, 
the components showed the following proper motion: for 
MCI 0.08±0.12 mas yr" 1 , for MC2 1.22±0.24 mas yr~\ 
for MC3 0.34 ± 0.10 mas yr" 1 , and for MC4 0.25 ± 0.05 
mas yr -1 , which correspond to 0.17±0.25c, 2.56±0.50c, 
0.71 ± 0.21c, and 0.53 ± 0.11c, respectively. MCI, MC3, 
and MC4 showed the sub-luminal speeds, which are con- 
sistent with the results shown in previous works (e.g., 
Piner et al. 2010). In these five epochs, MC2 showed 
super-luminal motion. Recently, Lico et al. (2012) also 
performed VLB A observations of Mrk 421 every month 
in 2011. In their 15 GHz observations, four jet compo- 
nents are also detected, and their angular distance from 
the core was consistent with MCI, MC2, MC3, and MC4. 

7. DISCUSSION 



4 



Niinuma ct al. 



To examine the behaviour inside the sub-pc region of 
Mrk 421, we conducted multi epoch VLBI observations 
anticipating a super-luminal component that is ejected 
from the core region after a large flare. Instead of detect- 
ing a newly emerging component on the VLBI images, 
we detected the apparent super-luminal inward motion of 
the jet component JC1. Such motions of jet components 
have been reported by previous MOJAVE project studies 
in which Lister et al. (2009b) discussed the detection 
probability of super-luminal inward motions among all 
observations to be approximately 2%. We consider our 
detection of the inward motion of JC1 in Mrk 421 is sim- 
ilar to the one rarely seen in the MOJAVE samples. 

There is also a possibility that the positional fluctua- 
tions of the jet component close to the core showed its 
unusual motion in short span(Edwards & Piner 2002), 
however, the motion of JC1 showed systematic inward 
motion at a super-luminal speed rather than its posi- 
tional fluctuations. Therefore, here we discuss several 
possible ideas suggested in Kellermann et al. (2004) for 
explaining an apparent super-luminal inward motion of 
JC1. Among them, it is clear that the free-free absorp- 
tion effect of disk, torus, and ambient cold plasma (the 
case discussed in Kameno et al. 2001) seems less likely. 
Another possibility they indicated is a highly curved jet. 
Although such extreme bending has been reported in a 
small number of sources (see, e.g., Savolainen et al. 2006), 
previous higher resolution studies for Mrk 421 do indeed 
not show any evidence for extreme bending within a par- 
sec of the jet. And, the possibility of the bobbling of the 
core component with a short time span is also consider- 
able. However, it seems that JC1 moved in closer to the 
core systematically through the first to the last epochs. 
Therefore, below, we focus on the following scenario pro- 
posed by Kellerman: a newly born component (JC0) on 
16th Feb 2010, when the largest X-ray flare occurred, 
was in the core but it was not resolved by our beam 
and it made the centroid of the core shift toward the jet 
motion. Also, as seen In Figure 5, the motions of JC1 
and MC3 seem discontinuous between the third and the 
forth epoch of JVN observations (Apr 25 and May 12, 
respectively), and the third epoch of MOJAVE obser- 
vation (Jul 12). It is, however, possible to explain this 
motion as JC1 (or MC3) got back on the track confirmed 
by previous studies, because the newborn component as- 
sociated with the giant flare, which shifted the centroid 
of the core faded within approximately 0.2 years. There- 
fore, basically, we have considered that the motion and 
the position of JC1 detected by our observations were 
consistent with previous studies. 

Lastly, we discuss the Dopplcr factor and the Lorcntz 
factor of JC0 under the assumption mentioned above. 
The apparent velocity of v\ a = — 3.48 ± 0.97c de- 
rived from the linear fit to the positions of JC1 rela- 
tive to the core, is the relative velocity. Therefore as- 
suming that JC1 did not move or its apparent veloc- 
ity is sufficiently slow, we can consider that the JC0 
moved toward JC1 with the apparent velocity of v$ a ~ 
3.48 ± 0.97c. In Figure 6, we show the bulk Lorentz 
factor (r ) and the Doppler factor (So) of JC0 calcu- 
lated by setting a range of 2.51c < vo a < 4.45c, and 
1° < < 5°, respectively. Here we use the equations 
of T = l/v/I - £2,<5 = 1/[F (1 - /3 o cos0)], (where, 



A) = V /c = (3oa/[ y/l + /3g a Sin(6> + (/>)], fioa = «0a/c, (j) = 

cos _1 (l/-^/l + /?q q )). As the result of these calculations, 
S a = 11.03+^, and T = 6.11±g;|| are obtained for the 
median value of 9 = 3°. 

We consider that these values are reasonable, because 
our result does not require extreme values for the jet 
parameters and the Doppler factor Sq is almost consis- 
tent with the typical value of ~ 10 derived from the 
broad band (optical to TeV 7-ray) SED fit (eg. Kino 
et al. 2002; BLazejowski et al. 2005; Donnarumma ct 
al. 2009; Abdo et al. 2011). 

8. CONCLUSION 

For giant X-ray flare of Mrk 421 occurred on 2010 
February 16, we carried out the multi-epoch VLBI obser- 
vations at intervals of ~ 20 days using JVN array soon 
after the giant flare. Although no new-born component 
can be seen directly in the total intensity images ob- 
tained by our quick follow-up observations, the jet com- 
ponent JC1 located at ~ 1 mas north-west from the core 
showed its apparent inward motion at an apparent speed 
of 3.48 ± 0.97c. As the most plausible explanation, we 
discuss that the inward motion of JC1 at a super-luminal 
speed was possibly caused by the ejection of a new com- 
ponent JC0, which could not be resolved with our ob- 
servations. In this case, the obtained Doppler factor of 
JC0 is around 10 to 20, which is consistent with the ones 
typically estimated by model fittings of spectral energy 
distribution for this source. 

To confirm our result, we must carry out a quick follow- 
up or sufficient densely monitoring for such a large flare 
with a phase-referencing VLBI technique. By doing this, 
there is a possibility of identifying the actual motion of 
all components including the core region after the large 
flare. Also, a polarization observation such as Pincr & 
Edwards (2005), or Marscher et al. (2008) is also an 
effective means of identifying whether a VLBI component 
is produced by the flare, and to investigate its motion 
right after the flare occurs. 
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TABLE 1 

Large X-ray and/or VHE 7-ray flares of Mrk 421 and VLBI follow-up 



X-ray or VHE 7-ray flare Astrophysical VLBI observations 



Epoch 


pa 
peak.VHE 

[Crab] 


F b 

pcak.X 

[mCrab] 


Ref. 


First epoch 


Intervals 


Ref. 


1996 May 


~ 10 


< 30 


1, 2 


1996 May 




12 


2001 Jan-Mar 


~ 13 




3 


2001 Mar 


3-5 months 


12, 13 


2004 Mar-Apr 


~ 3 


~ 135 


4 








2006 Jun 




~ 92 


5, 6 


2006 Jun 




12 


2008 Mar 




~ 110 


7 


2008 May 


~ 2 months 


12 


2008 Jun 




~ 125 


8 


2008 Jun 


~ 2 months 


12 


2009 Nov 




~ 96 


9 


2009 Dec 


3 months 


12 


2010 Jan 




120±10 


10 


2010 Jan 


3 months 


12 


2010 Feb 


> 10 


164±17 


10, 11 


2010 Mar 


~ 20 days 


14 



a The peak flux at VHE 7-ray band 
b The peak flux at X-ray band 

Reference - 1: Gaidos et al. (1996), 2:Schubnell (1997), 3: Krennrich et al. (2002), 4: BLazejowski et al. (2005), 5: Lichti ct al. 
(2008), 6: Lichti ct al. (2006), 7: Krimm et al. (2008), 8: Donnarumma et al. (2009), 9: D'Ammando et al. (2009), 10: Isobe et 
al. (2010), 11: Ong (2010), 12: MOJAVE program (http://www.physics.purdue.edu/astro/MOJAVE/index.html), 13: Piner & Edwards 
(2005), 14: This study 
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Niinuma ct al. 




07 Mar 10 (201 0.1 8) 




1 mas (0.64 pc) 



01 Apr 10 (201 0.25) 





25 Apr 1 (201 0.31 ) 1 2 May 1 (201 0.36) 




E< ' 

Fig. 1. — VLBI images of Mrk 421 with natural weighting at 22 GHz. Images are arranged from top left to bottom right in order of epoch. 
The epoch is indicated on the top of each panel as "DD MMM YY (YYYY.YY)". Also, the scale of 1 mas for our images is indicated on 
the panel of 1st epoch, two arrows on the image of 1st epoch show the core and the jet feature. These images are convolved with a circular 
beam size (full width at half maximum: FWHM) of 0.8 mas (shown in the bottom-left corner), which corresponds to the typical minor axis 
of synthesized beam in our observations. Contours begin at 3 mjy beam -1 and increase in -1, and 2 n steps. The peak flux and the image 
rms (-Fpeafe mjy beam -1 , F rm s mjy beam -1 ) at each epoch are (257.6, 0.6), (253.8, 1.5), (236.8, 0.5), (278.9, 0.8), respectively. 




„ 1.8 
to 

CO 

E 

o 
o 

® 1.4 
E 

S 1.2 



CO 



CO 

< 



0.6 



JC1 



0.1 0.2 0.3 0.4 

Elapsed year from 201 0.0 



0.5 



Fig. 2. — JC1 positions and weighted linear fit to JC1 positions (dashed line). The position uncertainties of each data point are uniformly 
rescaled so that the reduced x 2 to be unity. As the result of least square fit, we obtained the proper motion of —1.66 ± 0.46 mas yr -1 . 
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2010.18 2010.25 2010.31 2010.36 




Fig. 3. — VLBI images of Mrk 421 produced by restoring model- fit components. The epoch and the scale of 1 mas for our images are 
indicated on the top of each image as year "YYYY.YY", and in the bottom of the panel of the 1st epoch, respectively. The interval of 
each image is proportional to the date when each observation was made. The filled squares represent the position of JCl relative to the 
core, and the solid and the dashed lines show the trajectory of the core and of JCl, respectively. The decrease in the separation between 
the core and JCl is shown by these two lines. The beam size of this figure is the same as Figure 1, and is shown in the bottom-left corner 
of the panel of the 1st epoch. Contours begin at 3 mjy beam -1 and increase in -1, and 2™ steps. 




11 Feb, 10(2010.11) 




2 0-2-4 

Relative RA [Milli Arcsec] 

Fig. 4. — VLBA image of Mrk 421 with natural weighting at 15 GHz (MOJAVE program). The epoch is indicated on the top of each 
panel as "DD MMM, YY (YYYY.YY)", and the filled squares represent the position of each component relative to the core. The beam 
size of this figure is the same as Figure 1, and is shown in the bottom-left corner of the panel of the 1st epoch. Contours begin at 0.9 mjy 
beam -1 and increase in 2 n steps. 



TABLE 2 

Model fit parameters at each epoch - JVN 22 GHz. 



Epoch 


Component 




Model parameter 




Distance from the core 






Maj 
[mas] 


Axial Ratio 


PA 
[deg] 


Flux Density 
[mJy] 


R 

[mas] 


e 

[deg] 


2010 Mar 07 


Core 
JCl 


0.23±0.02 
0.65±0.01 


< 0.27 
1.0 


-24.0±4.5 


263.7 
17.9 


1.28 ±0.05 


-37.6±1.5 


2010 Apr 01 


Core 
JCl 


29+ 003 
^i+ - 35 

U ' Oi -0.51 


< 0.35 
1.0 




261.5 
9.1 


1.26 ±0.17 


-60.4±7.5 


2010 Apr 25 


Core 
JCl 


0.18±0.02 
0.60±0.06 


n 07+0.15 
1.0 


-34.81™ 


242.2 
16.8 


1.14 ±0.04 


-38.3±1.4 


2010 May 12 


Core 
JCl 


0.26±0.02 
0.57±0.11 


< 0.27 
1.0 


-34.0±4.7 


288.5 
14.3 


0.97 ±0.06 


-41.5±2.7 
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Elapsed year from 201 0.0 

Fig. 5. — Angular distances of each component detected with our observations and MOJAVE, relative to the core. The dashed and 
the solid arrows indicate the date when the large X-ray flares occurred. Also, each solid line shows the results of the linear-fit to each 
component without JC1. 



Range of the Doppler factor and the bulk Lorentz factor 



OOCKXX, 

<xx>ococ^ 
<xxxxxxx-. 
-xxxxxxxx^ 



Doppler factor 6 . 
(j3 Ol =3.48±0.97) 

bulk Lorentz factor [~ 
()3 =3.48±0.97) 



IMf>6(>bi 




1 1.5 2 2.5 3 3.5 4 4.5 5 

Viewing angle 6 [deg] 

Fig. 6. — <5n and Ty derived from the apparent velocity of the newborn VLBI component JC0 with the viewing angle ranging 1° < < 5°. 
At the viewing angle of 9 = 1°, 5n and fn exceed 16.8, 8.6, respectively. 
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TABLE 3 

Model fit parameters at each epoch - VLBA 15 GHz (MO J AVE). 



Epoch 


Component 




Model paramet' 


sr 


Distance from the core 


Maj 
[mas] 


Axial Ratio 


PA 

[dcg] 


Flux Density 
[mjy] 


R 

[mas] 


e 

[dcg] 


2009 Dec 17 


Core 


0.10 




-zo.o 


225.8 








MC4 


0.24 


1.0 




23.3 


0.38 


-31.7 




MC3 


0.37 


1.0 


- 


16.0 


0.95 


-35.1 




MC2 


0.63 


1.0 


- 


8.3 


1.60 


-32.9 




MCI 


1.27 


1.0 


- 


9.5 


4.04 


-32.5 


2010 Feb 11 


Core 


0.10 


0.00 


-73.1 


279.2 








MC4 


0.24 


l.U 




40.9 


u.oy 






MC3 


0.42 


1.0 


_ 


14.9 


1.17 


-35.6 




MC2 


0.97 


1.0 


- 


10.5 


2.12 


-33.5 




MCI 


1.16 


1.0 


- 


8.0 


4.09 


-33.7 


2010 Jul 12 


Core 


0.11 


0.73 


-52.5 


250.9 








MC4 


0.38 


i n 

l.U 




35.4 


U.4D 


AO A 




MC3a 


0.16 


1.0 


_ 


11.3 


0.96 


-50.0 




MC3 


0.52 


1.0 


- 


23.0 


1.31 


-37.6 




MC2 


0.46 


1.0 


- 


3.6 


2.32 


-29.4 




MCI 


1.40 


1.0 


- 


9.9 


4.22 


-31.6 


2010 Oct 15 


Core 


0.14 


0.54 


-30.5 


265.2 


- 


- 




MC4 


0.29 


1.0 




22.7 


0.56 


-39.1 




MC3 


0.46 


1.0 




17.6 


1.26 


-37.6 




MC2 


1.15 


1.0 




8.9 


2.58 


-37.7 




MCI 


0.58 


1.0 




4.6 


4.00 


-29.6 


2010 Nov 29 


Core 


0.11 


0.73 


-51.9 


248.4 








MC4a 


0.27 


1.0 




56.6 


0.25 


-26.3 




MC4 


0.23 


1.0 




17.3 


0.67 


-39.5 




MC3 


0.59 


1.0 




16.3 


1.38 


-40.2 




MC2 


1.04 


1.0 




4.3 


3.08 


-42.9 




MCI 


1.13 


1.0 




9.6 


4.23 


-30.8 



